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Abstract 

The length dependence is computed of the resistance of a disordered normal- 
metal wire attached to a superconductor. The scaling of the transmission 
eigenvalue distribution with length is obtained exactly in the metallic limit, 
by a transformation onto the isobaric flow of a two-dimensional ideal fluid. 
The resistance has a minimum for lengths near l/T, with I the mean free path 
and r the transmittance of the superconductor interface. 
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The resistance of a disordered wire increases with its length. In the metallic regime 
(length L much less than the localization length) the increase is linear, up to relatively small 
quantum corrections. The linear scaling seems obviously true regardless of the nature of 
the contacts to the wire, which would just contribute an additive, L-independent contact 
resistance according to Ohm's law. This is correct if the contacts are in the normal state. 
The purpose of this paper is to demonstrate and explain the complete breakdown of the 
linear scaling if one of the contacts is in the superconducting state. The resistance Rns{L) 
of a disordered normal-metal (N) wire (with mean free path I) connected to a superconductor 
(S) by a tunnel barrier (with transmission probability T) has a minimum when L ~ l/Y. 

The existence of a resistance minimum, and its destruction by a voltage or magnetic 
field, was first proposed by Van Wees et al. [[J, to explain the sharp dip in the differential 
resistance discovered by Kastalsky et al. 0. A similar anomaly has since been observed in 
a variety of semiconductor-superconductor junctions |§, and also in computer simulations 
A non- monotonous L-dependence of i?NS is implicit in the work of Volkov, Zaitsev, 
and Klapwijk ||, while Hekking and Nazarov |7] obtained a monotonous -Rns(^) f° r a high 
tunnel barrier. Very recently, Nazarov ||, using a Green's function technique, has confirmed 
the results of Ref . , in a highly interesting paper which has some overlap with the present 
work. 

Our analysis builds on two theoretical results. The first result we use is a Landauer-type 
formula for the conductance Gns = 1/-Rns> 

4e 2 N / T \ 2 

^-E(^). <U 

The numbers T n e [0, 1] are the eigenvalues of the matrix product tt\ with t the N x N 
transmission matrix of the disordered normal wire plus tunnel barrier (N being the number 
of transverse modes at the Fermi energy Ep). Eq. (|T]) holds in the zero-temperature, zero- 
voltage, and zero-magnetic field limit. Terms of order (A/Ep) 2 are neglected (with A the 
superconducting energy gap), as well as contributions from disorder in the superconductor. 
[Disorder in the superconductor will increase the effective length L of the disordered region 
by an amount of the order of the superconducting coherence length.] 



The second result we use is a scaling equation [10 



xln|A-A'|, (2) 

where s = L/l. We have made the conventional change of variables from T n to X n = 
(1 — T n )/T n , with \ n e [0, oo). The density p of the A's is defined by p(X, L) = (J2 n S(X — A n ), 
with (■ ■ •) the ensemble average. Eq. (§) describes the scaling of the eigenvalue density with 
the length L of the disordered normal region. In this formulation the tunnel barrier at the 
NS interface appears as an initial condition, 

P (A,o) = iV5(A-(i-r)/r), (3) 

where T is the transmission probability of the barrier. (For simplicity we assume a mode- 
independent transmission probability, i.e. all T n 's equal to T when L = 0.) Once p is known, 
one can compute the ensemble averaged conductance (Gns) from Eq. ([[]), 



2 



Ap 2 roc 

(G NS ) = — / d\p(X,s)(l + 2X)- 2 . (4) 
n Jo 



The non-linear diffusion equation (0) was derived by Mello and Pichard [0 from a 



Fokker-Planck equation |TT|JT^| for the joint distribution function of all N eigenvalues, by 



integrating out N — 1 eigenvalues and taking the large- N limit. This limit restricts its 
validity to the metallic regime (N ^> L/l), and is sufficient to determine the leading order 
contribution to the average conductance, which is O(N). The weak-localization correction, 
which is 0(1), is neglected (T^j. These 0(1) corrections depend on the magnetic field B, 
whereas the O(N) contributions to p described by Eq. (^) are insensitive to B. However, 
the relationship (f|) between Gns and p holds only for B = 0. (For B ^ 0, Gns depends 
not only on the eigenvalues of tt\ but also on the eigenvectors ||.) A priori, Eq. @j holds 
only for a wire geometry (length L much greater than width W), because the Fokker-Planck 
equation from which it is derived fl2| requires L 3> W. Numerical simulations [14| indicate 



that the geometry dependence only appears in the 0(1) corrections, and that the O(N) 
contributions are essentially the same for a wire, square, or cube. 

Our method of solution is a variation on Carleman's method \TE\. We introduce an 
auxiliary function 

F(z,s) = fdA'^M (5) 
Jo z — A 

which is analytic in the complex 2-plane cut by the positive real axis. Furthermore, 

lim F(z, s) = N/z. (6) 

\z\^oo 

The function F has a discontinuity for z = A±ie (with A > and e a positive infinitesimal). 
The limiting values F±(\, s) = F(X ± ie, s) are 

7r d r°° 

F ± = ±-p(X,s) + — dX'p(X',s)ln\X-X'\. (7) 
i aA Jo 

Combination of Eqs. (0) and (0) gives 

N§- g (F + - F.) = -|^A(1 + X)(Fl - Fl), (8) 
which implies that the function 

T(z, s) = N^-F(z, s) + £-z(l + z)F 2 (z, s) (9) 

is analytic in the whole complex plane, including the real axis. Moreover, T —>■ for 
| z | — > oo, in view of Eq. @. We conclude that T = 0, since the only analytic function 
which vanishes at infinity is identically zero. 

It is convenient to make the mapping z = sinh £ of the z-plane onto the strip S in the 
(-plane between the lines y = and y = — tt/2, where ( = x + iy. The mapping is conformal 
if we cut the z-plane by the two halflines A > and A < — 1 on the real axis. On S we define 
the auxiliary function U = U x + iU y by 
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The equation T = now takes the form 

^-U(C,s) + U((,s)^U((,s) = 0, (11) 

which we recognize as Euler's equation of hydrodynamics: (U x ,U y ) is the velocity field in 
the (x, y) plane of a two-dimensional ideal fluid at constant pressure. 

Euler's equation is easily solved. For initial condition [/(£, 0) = t/o(C) the solution to 
Eq. flnp is 

U(C,s) = U ((-sU((,s)). (12) 

To restrict the flow to S we demand that both ( and £ — sU((, s) lie in 5. From [/ we obtain 
the eigenvalue density, first in the x-variables 

p{x, s) = (2N/7c)U y (x - ie, s), (13) 

and then in the A- variables (A = sinh 2 x), 

p(A, s) = p(x, s)\dx/d\\ = p(x, s)\ sinh2x| _1 . (14) 

Eqs. (|T2|)-(|T4l) represent the exact solution of the non-linear diffusion equation (Q), for 
arbitrary initial conditions. 

The initial condition (|3|) corresponds to 

U (0 = isinh2C(cosh 2 C-r- 1 )- 1 . (15) 

The solution of the implicit equation (|T2|) is plotted in Fig. [1] (solid curves), for F — 0.1 and 
several values of s = L/l. For s>l and x <C s it simplifies to 

x = |arccoshr — |rs(r 2 — l) 1 / 2 cos a, (16a) 
a = nsN^ 1 p(x, s), r = a(Ts sin cr) -1 , (16b) 

shown dashed in Fig. [IJ. Eq. (^) has been obtained independently by Nazarov M . For s = 
(no disorder), p is a delta function at x , where T = 1/ cosh 2 a;o. On adding disorder the 
eigenvalue density rapidly spreads along the x-axis (curve a), such that p < N/s for s > 0. 
The sharp edges of the density profile, so uncharacteristic for a diffusion profile, reveal the 
hydrodynamic nature of the scaling equation (§). The upper edge is at 



s + ±ln(s/T) + e>(l). (17) 



Since L/x has the physical significance of a localization length |12[], this upper edge corre- 
sponds to a minimum localization length £ min = L/x mauX of order /. The lower edge at x ni \ n 
propagates from xo to in a "time" s c = (1 — T)/T. For 1 C s < s c one has 

aw = |arccosh (s c /s) - |[1 - (s/s c ) 2 ] 1/2 . (18) 
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It follows that the maximum localization length £ max = L/x m i n increases if disorder is 
added to a tunnel junction. This paradoxical result, that disorder enhances transmission, 
becomes intuitively obvious from the hydrodynamic correspondence, which implies that 
p(x, s) spreads both to larger and smaller x as the fictitious time s progresses. When s = s c 
the diffusion profile hits the boundary at x = (curve c), so that x m i n = 0. This implies 
that for s > s c there exist scattering states (eigenfunctions of ttf) which tunnel through 
the barrier with near-unit transmission probability, even if T <C 1. The number N open of 
transmission eigenvalues close to one (socalled "open channels" [|17]]) is of the order of the 
number of x n 's in the range to 1 (since T n = 1/ cosh 2 x n vanishes exponentially if x n > 1). 
For s ^> s c (curve e) we estimate 

iV O pen^p(0,s)=iV( S + r- 1 )- 1 , (19) 

where we have used Eq. (|T6|). 

To test these analytical results we have carried out numerical simulations similar to 
those reported in Ref. JJ. The disordered normal region was modeled by a tight-binding 
Hamiltonian on a square lattice (lattice constant a), with a random impurity potential at 
each site (uniformly distributed between ±^Up>). The tunnel barrier was introduced by 
assigning a non-random potential energy Ub = 2.3 E-p to a single row of sites at one end 
of the lattice, corresponding to a mode-averaged transmission probability V = 0.18. The 
Fermi energy was chosen at Ep = 1.5 u from the band bottom (with uq = h 2 /2ma 2 ). We 
chose Up, between and 1.5 u , corresponding to s between and 11.7 fL6| . Two geometries 
were considered: L = W = 285a (corresponding to N = 119), and L = 285a, W = 75a 
(corresponding to N = 31). In Fig. || we compare the integrated eigenvalue density v{x, s) = 
N^ 1 Jq dx' p(x', s), which is the quantity following directly from the simulation. The points 
are raw data from a single sample. (Sample-to-sample fluctuations are small, because the 
self-averaging quantities |T2|].) The simulation unambiguously demonstrates the 
appearance of open channels ( 0) on adding disorder to a tunnel barrier, and is 

in good agreement with our analytical result (TT2]), without any adjustable parameters. No 
significant geometry dependence was found, as anticipated. 

The average resistance of the NS junction is obtained directly from the complex velocity 
field U on the imaginary axis. From Eqs. (|4]) and ( |T0D we find 

i? NS = (a/2iVe 2 )lim ( ^_ W4 (dU/dCy 1 

= (h/2Ne 2 )(s + Q- 1 ), (20a) 

Q = f =-^-r(l + sin0) - l) , (20b) 

S COS \ 1 s cos <p J 

where G (0, 7r/2) is determined by 

<f>[\ - |r(l -sin0)] = rscos0. (20c) 



For r C 1 (or s ^> 1) Eq. (|20cj ) simplifies to = Tscos^, hence Q = Tsin^, in precise 
agreement with Ref. 0]. The scaling of the resistance with length is plotted in Fig. |3|. For 
T = 1 the resistance increases monotonically with L. The ballistic limit L — ► equals 
h/4Ne 2 , half the Sharvin resistance of a normal junction because of Andreev reflection [JTSJ . 
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For T < 0.5 a resistance minimum develops, somewhat below L = l/Y. The resistance 
minimum is associated with the crossover from a quadratic to a linear dependence of i?NS 
on l/T. The two asymptotic dependencies are (for T 1 and s 3> 1): 

R m = (h/2Ne 2 )s- l Y- 2 , if Ys < 1, (21a) 
R m = {h/2Ne 2 ){s + T- 1 ), if rs»l, (21b) 

to be contrasted with the classical series resistance ||, 

i?^ ss = (h/2Ne 2 ){s + 2r~ 2 ), (22) 

which holds if phase coherence is destroyed by a voltage or magnetic field. Eq. (^) would 
follow from a naive application of Ohm's law to the NS junction, with the tunnel barrier 
contributing an additive, disorder-independent amount to the total resistance. The quadratic 
dependence on 1 /Y in Eqs. ( |21a|) and (|22|) is as expected for tunneling into a superconductor, 
being a two-particle process 0. The linear dependence on l/T in Eq. ( [21 b|) was first noted 
in numerical simulations 0, as a manifestation of "reflectionless tunneling": It is as if one 
of the two quasiparticles can tunnel into the superconductor without reflection. Comparison 
of Figs. [I] and |^ provides the explanation. The resistance minimum occurs when the lower 
edge of the density profile reaches x = (curve c in Fig. [I]), and signals the appearance 
of scattering states which can tunnel through the barrier with probability close to one. 
For Ys ^> 1 (curve e), i?NS is dominated by the N open transmission eigenvalues close to one. 
From Eqs. (|IJ) and ( |T9"D we estimate i?NS — h/e 2 N open = (h/Ne^^s + Y' 1 ), up to a numerical 
prefactor, consistent with the asymptotic result fl21b ). 

It is essential for the occurrence of a resistance minimum that 1/-Rns depends non- 
linearly on the transmission eigenvalues. Indeed, if we compute the normal-state resistance 
l/i?N = (2e 2 /h) J2 n T n from the eigenvalue density fll^D, we find the linear scaling = 
(h/2Ne 2 )(s + r _1 ) for all Y and s. The cross-over to a quadratic dependence on l/T can not 
occur in this case, because of the linear relation between and T n in the normal state. 

In summary, we have presented a scaling theory for the resistance of a normal- 
superconductor microbridge. The scaling of the density p of transmission eigenvalues with 
length L is governed by Euler's equation for the isobaric flow of a two-dimensional ideal fluid: 
L corresponds to time and p to the ^/-component of the velocity field on the x-axis, with 
L/x corresponding to the localization length. This hydrodynamic correspondence provides 
an explanation for the resistance minimum which is both exact and intuitive. 

We thank Yu. V. Nazarov for valuable discussions and R. A. Jalabert for advise on 
the numerical simulations. This work was supported by the Dutch Science Foundation 
NWO/FOM and by the European Community. 
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FIGURES 



FIG. 1. Eigenvalue density p(x, s) as a function of x (in units of s = L/l) for T = 0.1. Curves 
a,b,c,d,e are for s = 2,4,9,30,100, respectively. The solid curves are from Eq. (12), the dashed 
curves from Eq. ([!(]) . The resistance minimum is associated with the collision of the density profile 
with the boundary at x = 0, for s = s c = (1 — r)/r. 



FIG. 2. Comparison between theory and simulation of the integrated eigenvalue density 
u(x,s) = iV _1 Jq dx' p(x' , s), for T = 0.18. The labels a,b,c indicate, respectively, s = 0,0.7,11.7. 
Solid curves are from Eq. (|l2|), data points are the x n 's from the simulation plotted in ascending 
order versus n/N = v (filled data points are for a square geometry, open points for an aspect ratio 
L/W = 3.8). The theoretical curve for s = is a step function at xq = 1.5 (not shown). The inset 
shows the full range of x, the main plot shows only the small-x region, to demonstrate the disor- 
der-induced opening of channels for tunneling through a barrier. (Note that, since T = l/cosh 2 x, 
x near zero corresponds to near-unit transmission.) 



FIG. 3. Dependence of the resistance i?NS on the length L of the disordered normal region 
(shaded in the inset), for different values of the transmittance T of the NS interface. Solid curves 
are computed from Eq. (20), for T = 1,0.8,0.6,0.4,0.1 from bottom to top. For T <C 1 the dashed 
curve is approached, in agreement with Ref. pj. 



